/*-------------------------<-- Start of Description -->-------------------------*\
| PURPOSE:  This macro completes a statistical analysis of many variables and    |
|           prints the output on one page.  Basically, it is a PROC UNIVARIATE   |
|           with an output that looks like PROC MEANS.                           |
|--------------------------<-- End of Description -->----------------------------|
|--------------------------------------------------------------------------------|
|------------------------<-- Start of Files Created -->--------------------------|
| SYNTAX:   %univ (indata   = _dataset_name_,                      <optional>    |
|                  outdata    = _output_dataset_,                  <optional>    |
|                  var    = _variable_list_,                       <optional>    |
|                  round  = _rounding_descriptor_,                 <optional>    |
|                  by     = _by_variable_,                         <optional>    |
|                  stats  = _statistics_to_print_)                 <optional>    |
\-------------------------<-- End of Files Created-->---------------------------*/
%macro univ(indata=, var=, outdata=contents, round=.01, by=, stats=_STD_) / des='Univariate statistics on one page';
%let notes=%sysfunc(getoption(notes, keyword));

options nonotes;
%let starttime=%sysfunc(datetime());
%local i dataset v round by bye pageby var fn stats bn lb;

/*----------------\
| assign defaults |
\----------------*/
%let the_end = N;

%if (%quote(&indata) = %quote()) %then %do;
   %let indata = %upcase(%scan(&sysdsn, 1).%scan(&sysdsn, 2));
   options &notes;
   %put NOTE: No datasets are specified.  The dataset &indata will be used.;
   options nonotes;
%end;
%else %let indata = %data(&indata);

%if (%quote(&var) = %quote()) %then %do;
   proc contents data=&indata out=contents(keep=name type) noprint; run;
   data _null_;
      set contents (where = (type=1));
      call symput ('var'||trim(left(_n_)), trim(left(name)));
      call symput ('fn', trim(left(_n_)));
   run;
   %do i = 1 %to &fn;
      %let var = &var &&var&i;
   %end;
%end;
%else %let fn = %wordcnt(&var);

/*------------------------------------------------------------\
| Calculate value of BY and PAGEBY--if by value is not %str() |
\------------------------------------------------------------*/
%if (%quote(&by) ne %quote()) %then %do;
   %let bn     = %words(&by);
   %let lb     = %scan(&by, &bn, %str( ));
   %let pageby = %str(pageby &lb ;);
   %let bye    = %str(by &by ;);
%end;

/*--------------------------------------------\
| If keywards are not used in the macro call, |
| this section assigns keywords based on      |
| option _STD_ (standard) or _ALL_ (all).     |
\--------------------------------------------*/
%if %upcase(&stats) = _STD_ %then %let stats = n mean std cv sum min p10 p25 p50 p75 p90 max;
%else %if %upcase(&stats) = _ALL_ %then
   %let stats = n nobs nmiss mean mode sum std cv css uss qrange
                msign normal stdmean probn probm probs signrank t
                kurtosis skewness range qrange sumwgt min p1 p5 p10
                p25 p50 p75 p90 p95 p99 max;

/*-------------------------------------\
| Use PROC CONTENTS to get label names |
\-------------------------------------*/
proc contents data=&indata(keep=&var &by) out=contents(keep=name label) noprint; run;

%let data_id = %sysfunc(open(&indata, i));
%do I = 1 %to &fn;
   %let label&i = %sysfunc(varlabel(&data_id, %sysfunc(varnum(&data_id, %scan(&var, &i, %str( ))))));
%end;
%let rc = %sysfunc(close(&data_id));

/*-------------------------------------------------\
| Univariate step writes out to dataset macro.temp |
\-------------------------------------------------*/
proc univariate data=&indata noprint; &bye var &var;
   output out      = temp
          n        = %_prefix(prefix=n,     from=1, to=&FN)
          mean     = %_prefix(prefix=mean,  from=1, to=&FN)
          std      = %_prefix(prefix=std,   from=1, to=&FN)
          cv       = %_prefix(prefix=cv,    from=1, to=&FN)
          min      = %_prefix(prefix=min,   from=1, to=&FN)
          p1       = %_prefix(prefix=p1_,   from=1, to=&FN)
          p5       = %_prefix(prefix=p5_,   from=1, to=&FN)
          p10      = %_prefix(prefix=p10_,  from=1, to=&FN)
          q1       = %_prefix(prefix=p25_,  from=1, to=&FN)
          median   = %_prefix(prefix=p50_,  from=1, to=&FN)
          q3       = %_prefix(prefix=p75_,  from=1, to=&FN)
          p90      = %_prefix(prefix=p90_,  from=1, to=&FN)
          p95      = %_prefix(prefix=p95_,  from=1, to=&FN)
          p99      = %_prefix(prefix=p99_,  from=1, to=&FN)
          max      = %_prefix(prefix=max,   from=1, to=&FN)
          sum      = %_prefix(prefix=sum,   from=1, to=&FN)
          kurtosis = %_prefix(prefix=kurt,  from=1, to=&FN)
          skewness = %_prefix(prefix=skew,  from=1, to=&FN)
          css      = %_prefix(prefix=css,   from=1, to=&FN)
          nmiss    = %_prefix(prefix=nmiss, from=1, to=&FN)
          msign    = %_prefix(prefix=msign, from=1, to=&FN)
          nobs     = %_prefix(prefix=nobs,  from=1, to=&FN)
          qrange   = %_prefix(prefix=qrang, from=1, to=&FN)
          t        = %_prefix(prefix=t,     from=1, to=&FN)
          normal   = %_prefix(prefix=norm,  from=1, to=&FN)
          sumwgt   = %_prefix(prefix=sumwt, from=1, to=&FN)
          stdmean  = %_prefix(prefix=stdmn, from=1, to=&FN)
          probn    = %_prefix(prefix=probn, from=1, to=&FN)
          range    = %_prefix(prefix=range, from=1, to=&FN)
          var      = %_prefix(prefix=var,   from=1, to=&FN)
          probm    = %_prefix(prefix=probm, from=1, to=&FN)
          mode     = %_prefix(prefix=mode,  from=1, to=&FN)
          signrank = %_prefix(prefix=sign,  from=1, to=&FN)
          uss      = %_prefix(prefix=uss,   from=1, to=&FN)
          probs    = %_prefix(prefix=probs, from=1, to=&FN) ;
run ;
/*--------------------------------------------------------------\
| transpose the dataset macro.temp to be printed, round values, |
| add labels, format the label descriptor, and keep variables   |
\--------------------------------------------------------------*/
data &outdata (keep = varname label min p1 p5 p10 p25 p50 p75 p90 p95
               p99 max kurtosis n skewness css nmiss std msign mean nobs
               sum qrange t normal sumwgt stdmean probn range cv probm
               mode signrank uss probs &by);
   set temp;
   label label    = 'Variable description'
         varname  = 'Variable name'
         n        = 'Number non-missing values'
         mean     = 'Mean value'
         sum      = 'Sum of all values'
         std      = 'Standard deviation'
         cv       = 'Coefficient of variation'
         min      = 'Minimum value'
         p1       = '1st %tile'
         p5       = '5th %tile'
         p10      = '10th %tile'
         p25      = '25th %tile'
         p50      = '50th %tile'
         p75      = '75th %tile'
         p90      = '90th %tile'
         p95      = '95th %tile'
         p99      = '99th %tile'
         max      = 'Maximum value'
         kurtosis = 'Kurtosis'
         skewness = 'Skewness'
         css      = 'Corrected sum of squares'
         nmiss    = 'Number missing values'
         msign    = 'Sign statistic'
         nobs     = 'Number observations'
         qrange   = 'Inter-quartile range'
         T        = "Student's T"
         normal   = 'Normality test statistic'
         sumwgt   = 'Sum of the weights'
         stdmean  = 'Standard mean'
         probn    = 'P-value of normality test statistic'
         range    = 'Range of values'
         var      = 'Variance'
         probm    = 'P-value of sign statistic'
         mode     = 'Most frequent value'
         signrank = 'Signed Rank Statistic'
         uss      = 'Uncorrected sum of squares'
         probs    = 'P-value of signed ranked statistic' ;
   length label $40 varname $8 ;
   /*------------------------------------\
   | This section mimmics PROC TRANSPOSE |
   \------------------------------------*/
   %do I = 1 %to &fn;
      varname  = "%scan(&var, &i, %str( ))";
      label    = "&&label&i";
      if label = '' then label = left(varname);
      n        = n&i;
      mean     = round(mean&i,  &round);
      sum      = round(sum&i,   &round);
      std      = round(std&i,   &round);
      cv       = round(cv&i,    &round);
      min      = round(min&i,   &round);
      P1       = round(P1_&i,   &round);
      P5       = round(P5_&i,   &round);
      P10      = round(P10_&i,  &round);
      P25      = round(P25_&i,  &round);
      P50      = round(P50_&i,  &round);
      P75      = round(P75_&i,  &round);
      P90      = round(P90_&i,  &round);
      P95      = round(P95_&i,  &round);
      P99      = round(P99_&i,  &round);
      max      = round(max&i,   &round);
      kurtosis = round(kurt&i,  &round);
      skewness = round(skew&i,  &round);
      css      = round(css&i,   &round);
      nmiss    = round(nmiss&i, &round);
      msign    = round(msign&i, &round);
      nobs     = round(nobs&i,  &round);
      qrange   = round(qrang&i, &round);
      T        = round(T&i,     &round);
      normal   = round(norm&i,  &round);
      sumwgt   = round(sumwt&i, &round);
      stdmean  = round(stdmn&i, &round);
      probn    = round(probn&i, &round);
      range    = round(range&i, &round);
      var      = round(var&i,   &round);
      probm    = round(probm&i, &round);
      mode     = round(mode&i,  &round);
      signrank = round(sign&i,  &round);
      uss      = round(uss&i,   &round);
      probs    = round(probs&i, &round);
      output ;
   %end;
run;

/*------------------------------------------------------\
| If by-variables are used, sort the transposed dataset |
\------------------------------------------------------*/
%if &by ne %str() %then %do ;
   proc sort data=&outdata; by &by; run;
%end ;

/*-----------------------------------\
| PROC PRINT is used to print output |
\-----------------------------------*/
proc print data=&outdata noobs label; &bye &pageby id varname; var label &stats; run;

/*-----------------------------------------\
| Delete temporary datasets and data views |
\-----------------------------------------*/
proc datasets lib=macro nolist memtype=(data view); delete contents view1 temp; run; quit;

/*--------------------------------------\
| Prepare the SAS session for open code |
\--------------------------------------*/
%put NOTE: Macro UNIV printed statistics for variabes %upcase(&var).;

%if (%quote(%upcase(&outdata)) ne %quote(CONTENTS)) %then %do ;
   %let outobs = %nobs(&outdata);
   %let outvar = %nvars(&outdata);
   options &notes;
   %put NOTE: The dataset %data(&outdata) has &outobs observations and &outvar variables.;
   options nonotes;
%end ;
%timenote (macro=univ,starttime=&starttime)
options &notes;
%mend univ;
<!--'"</title></head>-->

<script type="text/javascript">
//OwnerIQ
var __oiq_pct = 50;
if( __oiq_pct>=100 || Math.floor(Math.random()*100/(100-__oiq_pct)) > 0 ) {
var _oiqq = _oiqq || [];
_oiqq.push(['oiq_addPageBrand','Lycos']);
_oiqq.push(['oiq_addPageCat','Internet > Websites']);
_oiqq.push(['oiq_addPageLifecycle','Intend']);
_oiqq.push(['oiq_doTag']);
(function() {
var oiq = document.createElement('script'); oiq.type = 'text/javascript'; oiq.async = true;
oiq.src = document.location.protocol + '//px.owneriq.net/stas/s/lycosn.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(oiq, s);
})();
}
//Google Analytics
var _gaq = _gaq || [];
_gaq.push(['_setAccount','UA-21402695-19']);
_gaq.push(['_setDomainName','tripod.com']);
_gaq.push(['_setCustomVar',1,'member_name','duozhou',3]);
_gaq.push(['_trackPageview']);
(function() {
var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true;
ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s);
})();
//Lycos Init
function getReferrer() {
var all= this.document.cookie;
if (all== '') return false;
var cookie_name = 'REFERRER=';
var start = all.lastIndexOf(cookie_name);
if (start == -1) return false;
start += cookie_name.length;
var end = all.indexOf(';', start);
if (end == -1) end = all.length;
return all.substring(start, end);
}
function getQuery() {
var rfr = getReferrer();
if (rfr == '') return false;
var q = extractQuery(rfr, 'yahoo.com', 'p=');
if (q) return q;
q = extractQuery(rfr, '', 'q=');
return q ? q : "";
}
function extractQuery(full, site, q_param) {
var start = full.lastIndexOf(site);
if (start == -1) return false;
start = full.lastIndexOf(q_param);
if (start == -1) return false;
start += q_param.length;
var end = full.indexOf('&', start);
if (end == -1) end = full.length;
return unescape(full.substring(start, end)).split(" ").join("+");
}
function generateHref(atag, template){
atag.href=template.replace('_MYURL_', window.location.href.replace('http://', '')).replace('_MYTITLE_','Check%20out%20this%20Tripod%20Member%20site!'); 
}
var lycos_ad = Array();
var lycos_onload_timer;
var cm_role = "live";
var cm_host = "tripod.lycos.com";
var cm_taxid = "/memberembedded";
var tripod_member_name = "duozhou";
var tripod_member_page = "duozhou/MyOwnSetting/sasmacro/univ.sas";
var tripod_ratings_hash = "1714841670:7a08f44b345d526322806494c35b2e7d";

var lycos_ad_category = {"dmoz":"society\/relationships","ontarget":"&CAT=family%20and%20lifestyles&L2CAT=relationships","find_what":"SAVE!"};

var lycos_ad_remote_addr = "209.202.244.9";
var lycos_ad_www_server = "www.tripod.lycos.com";
var lycos_ad_track_small = "";
var lycos_ad_track_served = "";
var lycos_search_query = getQuery();
</script>

<script type="text/javascript" src="https://scripts.lycos.com/catman/init.js"></script>

<script type='text/javascript'>
 var googletag = googletag || {};
 googletag.cmd = googletag.cmd || [];
 (function() {
   var gads = document.createElement('script');
   gads.async = true;
   gads.type = 'text/javascript';
   var useSSL = 'https:' == document.location.protocol;
   gads.src = (useSSL ? 'https:' : 'http:') +
     '//www.googletagservices.com/tag/js/gpt.js';
   var node = document.getElementsByTagName('script')[0];
   node.parentNode.insertBefore(gads, node);
 })();
</script>

<script type='text/javascript'>
 googletag.cmd.push(function() {
   googletag.defineSlot('/95963596/TRI_300X250_dfp', [300, 250], 'div-gpt-ad-1450204159126-0').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_above_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-1').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_below_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-2').addService(googletag.pubads());
   googletag.pubads().enableSingleRequest();
   googletag.enableServices();
 });
</script>


<script type="text/javascript"> 
(function(isV)
{
    if( !isV )
    {
        return;
    }
    var adMgr = new AdManager();
    var lycos_prod_set = adMgr.chooseProductSet();
    var slots = ["leaderboard", "leaderboard2", "toolbar_image", "toolbar_text", "smallbox", "top_promo", "footer2", "slider"];
    var adCat = this.lycos_ad_category;
    adMgr.setForcedParam('page', (adCat && adCat.dmoz) ? adCat.dmoz : 'member');
    if (this.lycos_search_query)
    {
        adMgr.setForcedParam("keyword", this.lycos_search_query);
    } 
    else if(adCat && adCat.find_what)
    {
        adMgr.setForcedParam('keyword', adCat.find_what);
    }
    
    for (var s in slots)
    {
        var slot = slots[s];
        if (adMgr.isSlotAvailable(slot))
        {
            this.lycos_ad[slot] = adMgr.getSlot(slot);
        }
    }

    adMgr.renderHeader();
    adMgr.renderFooter();
}((function() {

var w = 0, h = 0, minimumThreshold = 300;

if (top == self)
{
    return true;
}
if (typeof(window.innerWidth) == 'number' )
{
    w = window.innerWidth;
    h = window.innerHeight;
}
else if (document.documentElement && (document.documentElement.clientWidth || document.documentElement.clientHeight))
{
    w = document.documentElement.clientWidth;
    h = document.documentElement.clientHeight;
}
else if (document.body && (document.body.clientWidth || document.body.clientHeight))
{
    w = document.body.clientWidth;
    h = document.body.clientHeight;
}
return ((w > minimumThreshold) && (h > minimumThreshold));
}())));




window.onload = function()
{
    var f = document.getElementById("FooterAd");
    var b = document.getElementsByTagName("body")[0];
    b.appendChild(f);
    f.style.display = "block";
    document.getElementById('lycosFooterAdiFrame').src = '/adm/ad/footerAd.iframe.html';
    


    
    // DOM Inj Ad
    (function(isTrellix)
    {
        var e = document.createElement('iframe');
        e.style.border = '0';
        e.style.margin = 0;
        e.style.display = 'block';
        e.style.cssFloat = 'right';
        e.style.height = '254px';
        e.style.overflow = 'hidden';
        e.style.padding = 0;
        e.style.width = '300px';


        var isBlokedByDomain = function( href )
        {
            var blockedDomains = [
                "ananyaporn13000.tripod.com",
                "xxxpornxxx.tripod.com"
            ];
            var flag = false;
            
            for( var i=0; i<blockedDomains.length; i++ )
            {
                if( href.search( blockedDomains[ i ] ) >= 0 )
                {
                    flag = true;
                }
            }
            return flag;
        }

        var getMetaContent = function( metaName )
        {
            var metas = document.getElementsByTagName('meta');
            for (i=0; i<metas.length; i++)
            { 
                if( metas[i].getAttribute("name") == metaName )
                { 
                    return metas[i].getAttribute("content"); 
                } 
            }
            return false;
        }
        
        var getCommentNodes = function(regexPattern)
        {
            var nodes = {};
            var nodesA = [];
            var preferredNodesList = ['a', 'c', 'b'];
        
            (function getNodesThatHaveComments(n, pattern)
            {
                if (n.hasChildNodes())
                {
                    if (n.tagName === 'IFRAME')
                    {
                        return false;
                    }
                    for (var i = 0; i < n.childNodes.length; i++)
                    {
                        if ((n.childNodes[i].nodeType === 8) && (pattern.test(n.childNodes[i].nodeValue)))
                        {
                            var areaName = pattern.exec(n.childNodes[i].nodeValue)[1];
                            nodes[areaName] = n;
                        }
                        else if (n.childNodes[i].nodeType === 1)
                        {
                            getNodesThatHaveComments(n.childNodes[i], pattern);
                        }
                    }
                }
            }(document.body, regexPattern));

            for (var i in preferredNodesList)
            {
                if (nodes[preferredNodesList[i]])
                {
                    if( isTrellix && nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode )
                    {
                        nodesA.push(nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode);
                    }
                    else
                    {
                        nodesA.push( nodes[preferredNodesList[i]] );
                    }
                }
            }
            return nodesA;
        }
        
        
        var properNode = null;
        var areaNodes = getCommentNodes( new RegExp( '^area Type="area_(\\w+)"' ) );

        for (var i = 0; i < areaNodes.length; i++)
        {
            var a = parseInt(getComputedStyle(areaNodes[i]).width);
            if ((a >= 300) && (a <= 400))
            {
                properNode = areaNodes[i];
                break;
            }
        }


        var propertyName = getMetaContent("property") || false;
        if( isTrellix && (properNode) )
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            properNode.insertBefore(e, properNode.firstChild);
        }
        else if( isTrellix && !( properNode ) ) // Slap the ad eventhought there is no alocated slot
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            e.style.cssFloat = 'none';
            var cdiv = document.createElement('div');
            cdiv.style = "width:300px;margin:10px auto;";
            cdiv.appendChild( e );
            b.insertBefore(cdiv, b.lastChild);
        }
        else if( !isBlokedByDomain( location.href ) )
        {
            var injF = document.createElement('iframe');
            injF.style.border = '0';
            injF.style.margin = 0;
            injF.style.display = 'block';
            injF.style.cssFloat = 'none';
            injF.style.height = '254px';
            injF.style.overflow = 'hidden';
            injF.style.padding = 0;
            injF.style.width = '300px';
            injF.src = '/adm/ad/injectAd.iframe.html';

            if( b && ( !isTrellix || ( typeof isTrellix == "undefined" ) ) ) // All other tripod props
            {
                var cdiv = document.createElement('div');
                cdiv.style = "width:300px;margin:10px auto;";
                cdiv.appendChild( injF );
                b.insertBefore(cdiv, b.lastChild);
            } 
        }
  }( document.isTrellix ));
}

</script>

<div id="tb_container" style="background:#DFDCCF; border-bottom:1px solid #393939; position:relative; z-index:999999999!important">
<!--form name="search" onSubmit="return searchit()" id='header_search' >
<input type="text" placeholder="Search" size=30 name="search2" value="">
<input type="button" value="Go!" onClick="searchit()">
</form>
<style>
form#header_search {
    width: 916px;
    margin: 0 auto 8px;
    position: relative;
}


form#header_search input {
    height: 40px;
    font-size: 14px;
    line-height: 40px;
    padding: 0 8px;
    box-sizing: border-box;
    background: #F4F2E9;
    border: 1px solid #BBB8B8;
    transition: background-color 300ms ease-out,
                color 300ms ease;
}

form#header_search input[type="text"] {
    width: 100%;
}
form#header_search input[type="text"]:focus {
    border-color: #A2D054;
    background-color: #fff;
    box-shadow: 0 0px 12px -4px #A2D054;
}



form#header_search input[type="button"] {
    position: absolute;
    top: 1px;
    right: 1px;
    opacity: 1;
    background: #DFDCCF;
    color: #463734;
    width: 125px;
    cursor: pointer;
    height: 38px;
    border: none;
}
form#header_search input[type="text"]:focus ~ input[type='button']:hover,
form#header_search input[type='button']:hover {
    background-color: #A5CE56;
    color: #fff;
}
form#header_search input[type="text"]:focus ~ input[type='button'] {
    background-color: #52AEDF;
    color: #fff;
}

</style>

<script>
function searchit(){
    
    // determine environment 
    var search_env 
    if (lycos_ad_www_server.indexOf(".pd.") > -1) {
       search_env = 'http://search51.pd.lycos.com/a/';
    } else if (lycos_ad_www_server.indexOf(".qa.") > -1) {
       search_env = 'http://search51.qa.lycos.com/a/';
    } else {
       search_env = 'http://search51.lycos.com/a/';
    }

var search_term = encodeURIComponent(document.search.search2.value)
var search_url = search_env+search_term;
window.open(search_url);

return false
}
</script-->
<style>
    .adCenterClass{margin:0 auto}
</style>
<div id="tb_ad" class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">

<div id="ad_container" style="display:block!important; float:left; width:728px ">
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "leaderboard" in lycos_ad) {
  document.write(lycos_ad['leaderboard']);
}
</script>
</div>
</div>
</div>
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "slider" in lycos_ad) {
  document.write(lycos_ad['slider']);
}
</script> <!-- added 7/22 -->
<div id="FooterAd" style="background:#DFDCCF; border-top:1px solid #393939; clear:both; display:none; width:100%!important; position:relative; z-index:999999!important; height:90px!important"> 
<div class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">
<div id="footerAd_container" style="display:block!important; float:left; width:728px">
<iframe id="lycosFooterAdiFrame" style="border:0; display:block; float:left; height:96px; overflow:hidden; padding:0; width:750px"></iframe>
</div>
</div>
</div>